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TECHNICAL MEMORANDUM X-34J49 

TWO-DIMENSIONAL CONVOLUTE INTEGERS FOR 
OPTICAL IMAGE DATA PROCESSING AND SURFACE FITTING 

INTRODUCTION 


Two-dimensional convolute integers are sets of convolution coefficients 
of integer value which can be used in optical image data processing and surface 
fitting problems. Applying these coefficients to a set of data is equivalent to 
doing two-dimensional regression calculations; i.e. , fitting a surface by least 
squares, without considering the time-consuming summations associated with 
the usual normal equations. The only requirement for use of these coefficients 
is that the data points on the surface be either equidistant along their orthogonal 
axis, or be spaced an integer multiple of a common factor. 

A moving smoothing average is one of the fastest algorithms which can 
be applied for noise filtering in an optical image, but it tends to reduce reso- 
lution. The coefficients in such a filter are convolution coefficients. Two- 
dimensional convolute integers utilize the fast logic of the moving smoothing 
average without sacrificing resolution. From the definition of these coefficients, 
ten properties appear. These properties indicate the existence and universal 
aspect of the filter numbers. As a consequence of being derived from regres- 
sion theory, the filters have low-pass and high-pass filtering properties. This 
is equivalent to smoothing and generating the partial derivatives associated 
with surface order. The number of unique integers per filter represents the 
minimum number of multiplications per image point and significantly affects 
the filtering time. A zero order or first order surface is identical to nearest 
neighbor averaging. Double convolution, which represents band-pass filtering 
and can generate surface normals, can be accomplished by a single set of 
integers using an expanded filter mask. 

The appendices contain detailed calculations for a 23-point filter, 

5 by 5 mask, second order surface. The normal equations for a first and 
third order are also presented. 
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MOVING SMOOTHING AVERAGE 


The concept of a moving smoothing average represents one of the 
simplest and fastest filters that can be applied to noise reduction in a two- 
dimensional image. The filter can be considered a mask, m rows by n 
columns, which when applied to the data in a two-dimensional image causes a 
nearest neighbor interaction. 

The interaction results from adding the intensities of all the points sur- 
rounding a particular location, including the value of the intensity at this 
location in the sum, and then dividing by the total number of points considered. 
In this way, nearest neighbors have an opportunity’ to influence the new average 
value, which shall be considered the filter data point. 

This m by n mask can be passed over all the data points in an image, 
with the exception of a border of points left unfiltered. Requiring an odd 
number of rows and columns unambiguously determines the center point of the 
filter mask, which can then be shifted column by column, advancing one row at 
a time when a border column is reached. 


In particular, each data point in the mask can lx? thought of as multi- 
plied by a coefficient C lEq. (l) j, added with all other points in the mask, 

the sum divided by a normalizer (Norm), and the resulting new filtered data 
point put in place of the raw data point. By this method, the raw data array is 
transformed into a filtered data array. 


For a moving smoothing average, all the C., terms equal unity, and 


Norm is equal to the number of points in the array. 


z' 

rc 


(x,y) 


m n 

X X C z (x,y)/Norm 
i--m j -n 1 1 


( 1 ) 


where 

z - array of raw data points, 
z' = array of filter data points, 

ij represents a point in the (2m h 1) by (2n •< 1) filter, and rc 
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represents a raw data point in an image at row r and column c and is the 
center point of the filter mask. 

Such a filter is nonrecursive and, therefore, requires raw data points 
in the Z vector at all times. Unfortunately, this type of filtering, though simple 
and fast, tends to blur an image by decreasing the overall sharpness, or reso- 
lution. 


The problem which then arises requires the development of a set of 
weighting coefficients C _ and Norm which filter with the speed of Eq. ( 1) and 

yet do not decrease resolution. 


CONVOLUTION 

The coefficients C _ in Eq. ( 1) are convolution coefficients for the 

following reason. The filter can be considered an operator which forms the 
filtered data by integrating the raw data over a weighting function Q(o,p ): 

z'(x.y) = Jn(a,/3)z(x-a, y-p)dadp . (2) 

This integral is defined as the convolution of z(x-o; y-fi) with Q(a,p ). In a 
digital filter, the weighting function is of the form 


n(a,p) 


m 


Z 

i=-m 


n 

Z + i, P + j)/Norm 

j=-n J 


(3) 


where 6 is the Dirac delta function representing the discrete sampling of the 
data. Using this in Eq. (2) gives 


z'(x,y) 


/ 


m 

£ 

i=-m 


n 

V C 6( a + i, p + j)/Norm 
j=-n lJ 


z(x-a, y~p)dadp 


m n 

= Z Z C z( x + i, y + j)/Norm 

i=-m j=-n J 


(4) 
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When z(x+ i, y+ j) is equated with z (x,y), Eq. (4) is seen to be identical 

-J 

to Eq. (1). Thus, the C.^'s and Norm are known as convolution coefficients. 

TWO-DIMENSIONAL REGRESSION 


Two-dimensional regression calculations are equivalent to fitting a 
polynomial surface to a set of data by the method of least squares [l], A 
general expression for a two-dimensional surface is 


m 

- Z 

i=0 


n 


j-o 


A..x| yf 
ij k J k 


(5) 


Equation (5) represents the intensity z\ at the point k on the surface of 

K 

order m 4 n described by the regression coefficients A„. 

Let z^ represent the actual intensity at the kth position on an image. 
Fitting the surface to the data set z^(x,y) by the residual sum of squares 
leads to Eq. ((j), 


6 2 = 




Z 

j=0 


Vk yi k 


(fi) 


Selecting the best fit regression coefficients, 


np / m n \ 

a6VaA uv ■ 2 k z (\-L o l )vl^ 


(7) 


leads to the normal equation 


m n np np 

Z V v a l+u j+v V u v 

l L A ii x k \ = L \ \ y k 

i=0 j=0 k=l 1 k=l k k k 
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where u = 0,m and v = 0,n. Matrix representation of Eq. (8) results in the 
more familiar form of the normal equations. 

XI T • XY • A= XT T • Z (9) 

where Z is a column vector of raw data points, np by 1; A is a column vector 
of coefficients representing the surface, nt by 1; ST is a matrix of cross 
product terms, np by nt; and 


nt = (m + n+1) • (m4n+2)/2 


number of rows and columns in the symmetric XY 


T 


SY matrix. 


TWO-DIMENSIONAL CONVOLUTE INTEGERS 


Definition 

Rewriting ^q. (1) in matrix form allows an identification of the 
regression coefficients, A, with the convolution coefficients C. 


z' = C • Z/Norm (10) 

rc 


(Note that C is a row vector of double-subscripted elements, Z is a column 
vector of double subscripted elements, and both have np elements. ) 

Solving the normal equations for A, which is equivalent to obtaining intensity 
values on the fitted surface, yields 



Recognizing the similarity in form between Eq. (10) and ( 11) 

(the A's and Z's are both linear combinations of the Z's) permits a definition 
of the A's in terms of the C's, 


A 


U 


( 12) 


C • Z/Norm 
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Evaluation of the C's and Norm associated with each is readily 

accomplished by applying Cramer's Rule for obtaining the A's in Eq. (0) and 
identifying the coefficients of Z as C/Norm. 


Properties 

The convolution coefficients have ten properties which can be seen 
by solving Eq. (11) and identifying the C's and Norm. These properties, 
listed in Table 1, will be discussed in detail and lead to the application of 
the sets of C's and Norm's in optical image data processing and surface 
fitting. 


TABLE 1. PROPERTIES OF THE TWO-DIMENSIONAL 
CONVOLUTION COEFFICIENTS 

1. Existence - Number of points on surface must be 
greater than surface order. 

2. Universal Numbers ~ Convolution coefficients depend 
only on surface order and filter mask size. 

3. Antisymmetry The cross product matrix is sparse. 

4. Equal Interval — Convolution coefficients are integral 
numbers. 

5. Regression - Convolution coefficients derived from 
regression theory arc filters, 

6. Surface Fitting - Convolution coefficients represent 
surface partial derivatives. 

7. Equivalences - Convolution coefficients for adjacent 
surface orders are identical, as are the transpose 

of complementary roughing filters of the same order. 

8. Symmetries - Convolution coefficients have a minimum 
unique set of values per filter. 

9. Plane Surface — Convolution coefficients for a first 
order surface represent a moving smoothing average. 

10. Double Convolution - Convolution coefficients for 
band-pass filtering and surl. ■ e normals represent 
double convolution. 
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Existence . For each A associated with a surface of order m+n, there 
exist a set of C's and a value for Norm so long as 


SI • XI 4 0 


(13) 


This nonsingularity requirement indicates that the number of points on the 
surface, np, must lx; greater than the order of the surface, m+n. 

Universal Numbers . The values of C and the value of Norm associated 
with a surface of order nn n and np points are independent of the values of the 
intensities, z , on the surface. Thus, the convolution coefficients arc 

K 

universal numbers completely descrilx;d bv two numbers, m+n and np. 

Antisymmetry . An antisymmetric filter mask has an odd numlx*r of 
rows and columns, not necessarily equal, which, as previously stated, 
unambiguously establishes the center point as 0,0. Such a mask reduces the 
elements in the symmetric cross product matrix to zero where 

p q 

L \ y k ° * il4) 

k-1 k 


when either p or q is odd. Equation ( 14) is true since the single sum over all 
the data points in the filter mask can lx* replaced by a double summation. 

Eq. ( 1) or Eq. ( 15). 


np 

l 

k 1 


r 


Z 


c 

S' 

i-j 

y -c 


(15) 


where 


np = (2r+ 1) • 12c : 1) 

Substitution of this expression in Eq. ( 11) leads to the sparseness for anti- 
symmetric mask since each sum on the right-hand side of Eq. (10) is zeio. 
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r „ c 

2*2 y q 

=-r v=~f 


= 0 


( 16 ) 


when either p or q is odd. Because of the sparseness, the regression 
coefficients A„ can be partitioned into matrices of lower order (see Appendix A). 


Equal Interval . An equal interval filter mask has the spacing between 
data points an integer multiple of a common factor. The common factor between 
rows Ar need not equal the common factor between columns Ac, 


Ar^ Ac 


( 17 ) 


and the integers need not be consecutive but must remain fixed as the mask 
moves over an image (otherwise the C's and Norm change). As a consequence, 
the terms in the cross product matrix XT are integer values, exclusive of the 
common factors Ar and Ac, and accordingly the C's and Norm are integer 
values. 


Factoring common integers from the C's and Norm leads to a uni/ n rsal 
set of integers associated with each A., which can be dcscriltcd as two-dimen- 

sional convolute integers [2]. Note that the antisymmetry and equal interval 
properties are mutually independent. 


Regression . As a result of being derived from regression theory, sets 
of C's and Norm's satisfy two filtering principles. The A., terms can be 

obtained from the normal equations (see Appendix A) by Ciamcr's Rule. Now 
the coefficients of z on the right-hand side of the normal equations are identical 
to the first column of the XX * XX ^ matrix * sec Appendices A and B). Solving 
the normal equations for A., and factoring the Zz terms from the i+j column of 


the numerator clearly identifies the C's and Norm (see Appendix A). Summing 
the C's over all the data points in the mask and considering the determinant of 
the numerator leads to two filtering properties 1 3] , 


( 1) Smoothing: The set of C's and Norm associated with A satisfy 


np 

Yj C /Norm = 1 (is) 

k=l 
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low-pass filtering. 


(2) Roughing: The set of C's and Norm associated with A.., where 
i,i ^ 0, satisfy l -* 

np 

2 C /Norm = 0 ( 19 ) 

k-1 K 


band-pass or high-pass filtering. 

Surface Fitting . The two-dimensional convolution filters generate the 
various partial derivatives for a surface of order m+n fitted to a data set by 
regression. This can be seen by identifying all the regression coefficients in 
Eq. (5), evaluated at the center point of the mask. 

The constant term, A qo , is the fitted value of the intensity or a smoothed 

data point (the zero partial derivative). Thus, the C's and Norm associated 
with A represent smoothing filters. 

Expanding Eq. (5), assuming an antisyrr metric equal interval 
(Ar £ Ac) mask, taking partial derivatives evaluated nt the central point 
(x,y =0), leads to Eqs. (20) through (23). 






Ar Ac^ 9x*9y^ 


where 


( 20 ) 


r = x * Ar + a , 

row spacing; 

c = y • Ac + (j , 

column spacing; 


and x and y are antisymmetric integer values. At the center point in the filter 
mask, 





Thus, the C's and Norm's associated with A.. where i, j £ 0 represent partial 
derivative (roughing) filters. Note that Eq. (23) indicates that A terms 
are independent of the row column spacing ( Ar, Ac). 

Equivalences . As a result of the antisymmetric properties of the 
filter mask, a number of equivalences exist between convolution coefficients 
of adjacent surface orders and between the transpose of complementary 
roughing filters of the same surface order. The effect of these equivalences 
is to reduce the number of calculations required to present all the filters 
associated with the various surface orders and mask sizes. 

Table 2 lists the two types of equivalences in the convolution 

coefficients: C (ij), where £ represents surface order and ij the appropriate 
partial derivative. The proof of these properties is presented in Appendix B, 
which shows the normal equations for a first order and third order surface. 

Symmetries . Some general symmetry properties exist for the set of 
C's as dated with each A... These properties depend on ij odd or even and 

result from the symmetry associated with the normal equations, the anti- 
symmetry of the mask, and the equal interval between data points. The 
symmetry properties are expressed in Table 3. An alternate expression for 
the convolution coefficients leads to these properties, Appendix B. 


Coefficients which lie on a zero row or column have a two-fol- 
. symmetry. All other coefficients have at least a four-fold symmetry; and 

for a symmetric mask, nondiagonal, nonzero, row or column terms possess 
• eight-fold symmetry. 


s 

> 


i 

t 

t 






TABLE 2. EQUIVALENCES 


Regression 

Coefficient 


A ( ij ) 
i+j odd 
i+j even 


Complementary 

Transpose 


Surface Order 

f odd 

i even 

i f+i 

C = C 

i l-i 

C = C 

f i-i 

c = c 

f f+i 

c = c 

Identical Surface Order 

C(ji) * 

= c(ij) T 


a. f represents the order of the surface. 

b. ij represents the appropriate partial derivative. 

c. T represents transpose. 

























These properties allow the number of multiplications per filter to be 
co.' : iderably reduced. Table 4 indicates the number of multiplications per 
filter. In general this table gives the maximum number of multiplications, 
f.poufic filters are capable of having fewer unique numbers; e. g. , multipli- 
cations. These properties sufficiently affect the time required to filter an 
image. 

A Plane Surface . A satisfying, though possibly trivial result, is the 
equ- alence to be found between regression calculation filtering by a plane 
surl.ice and nearest neighbor averaging, which was discussed under the heading 
Moving Smoothing Average. Appendix B indicates the reduced form of the 
normal equations for a first order surface. For smoothing, the C*s are clearly 
unitj and Norm = np, the number of points considered. Smoothing an image by 
applying the two-dimensional convolute integers for a first order surface is 
identical to a moving smoothing average. 

Double Convolution . The convolute integers generated so far have 
satisfied low-pass and high-pass filtering criteria. Combining these filters 
v. ill result in band-pass filtering, or double convolution. 

Let 7.P represent the filtered value of a point in an image as a result of 
k 

coma lving the raw data with a pth order, rth row and cth column filter. 


?■, - ^ ) ) C P z /Norm 

k . 1 i Li Li mn mn t 

Ar Ac' m^-r n^-c 


This may be a smoothing (i,j = 0) or high-pass filter (i, j =£0). 

Assuming the data needs additional filtering and the mask spacing 

(Ar, A •) remains fixed, let represent convolving the already filtered data 

k 

with a surface of qth order. 


w . + n 

Ar* Ae* 


r' c 

Z Z 

m* — r* n' -c' 


c q , ,z q ,/Norm 
m 'n ' m 'n " q 
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TABLE 4. MULTIPLICATIONS PER FILTER ab 


Mask 

No. of Multiplications 

Symmetric, 


i = j 


even, even 

(np+1) (np+3)/8 

odd, odd 

(np-1) (np+l)/8 

Nonsymmetric, 


i * 3 


even, even 


i,j < order 

^(np-1) (np+l)/4j + 1 

i, j = order 
(j,i = 0) 

(np+l)/2 or (np-l)/2 

odd, odd 


i,j < order 

(np-l)/2 

i, j = order 0 


even, odd 


i , j < order 

[<np-l)/2] 2 or (np-1) (np+l)/4 

i, j = order 
(j.i= 0) 

(np-l)/2 


a. np represents the number of data points in a row or column 
of the filter mask. 

b. i,j represents the appropriate partial derivative. 

c. i, j combination not possible. 
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,qp _ (r + n , .(i + m 


y y y v r q 
i+i' A J+T p . L L L Sn'n‘ 

Ac * m'=-r' n'=-c m=-r n=-c 


♦ 2 . ./Norm Norm 

mn m+m n+n' p q 


Recognizing that the factorial terms and the powers of Ar and Ac are merely 
gain factors (which in no way affect the filtering) , the double convolution can be 
combined into a single convolution where 


rt-r (jtu 

■?- Z 2 

a=-(rfr“) p=-(c+c') 


C a z J Norm, 
at / 3 a$ 


where 


r c 


Jj # Jj „ Z Z C m'n- C i 

r m — - -h M — o 


m =-r n =-c m=-r n=-c 


and where the sums are constrained a - m'+ri and fi = n'+n, and where 


Norm = Norm Norm 

P q 


The C Q(i *s and Norm show that double convolution can be accomplished 

by a new single set of two-dimensional convolute integers using a[2(r+r') + l] 
by [2(c+c') + 1] expanded mask. A considerable savings in time results from 
applying the results of Eqs. (28) and (29). Indeed, double convolution 
represents the means by which a single filter can be constructed for obtaining 
the normal to the surface, at points in an image, or on a two-dimensional 
surface. 




CONC LUSION 


The theory of two-dimensional convolute integers has been discussed 
emphasizing the filtering aspects of the coefficients for optical image data 
processing. The coefficients generate the various partial derivatives 
associated with the surface order. Rand-pass filtering and surface normals 
can be obtained from a single set of integers with an expanded mask represen- 
ting double convolution. 

The coefficients have the filtering speed of a moving smoothing average 
and retain resolution since they represent two-dimensional regression calcu- 
lations. Being derived from regression theory, the coefficients satisfy the 
criteria for low-pass and high-pass filtering. 

Based upon a few simple concepts”” (l) antisymmetry, (2) equal 
interval, and (3) mask rows equal columns — , a number of equivalences exist 
between surface orders and among filter coefficients. All these properties tend 
to reduce die number of unique coefficients per filter; e.g. , the number of 
multiplications per filter point, or enhance the speed of the filter. The integer 
nature of the coefficients also tends to speed the filtering operation. Each 
filter is completely described by two numbers, mask size and surface order. 

Detailed calculations have shown the surface order and coefficient 
equivalences between first, second, and third order surfaces and indicate 
the manner by which two-dimensional convolute integers for any size mask 
and any order surface can be obtained. 

A companion report is planned which will deal with the software and 
hardware aspects of the filters and provide extensive tabulated data on filter 
coefficients. 
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TWO-DIMENSIONAL CONVOLUTE INTEGERS 
SECOND ORDER SURFACE 
5 BY 5 MASK 

The normal equations for the quadratic surface are 


00 

01 

02 

10 

11 

20 



uv 

np 

2y 

2y 2 

2 x 

2 xy 

^ ^ Aqq 


2 z 

00 

2y 

2y 2 

2y 3 

2 xy 

2xy 2 

2x2 y a 01 


2 zy 

01 

2y2 

Sy3 

Sy 4 

2 xy 2 

2 xy 3 

2x2 y 2 Agj 


2 zy 2 

02 

2x 

2xy 

2xy 2 

2 x 2 

2 x 2 y 

2 x 3 A io 


2 zx 

10 

Zxy 

Sxy 2 

2 xy 3 

2 x?y 

2xV 

2x 3 y A u 


2 zxy 

11 

Ex 2 

2x 2 y 

2xV 

2k 3 

2 x 3 y 

2 X 4 Ajq 


2 zx 2 

20 


For brevity, all summations are sums over all the raw data points k = 
l,np; and m + n is constrained to be £ 2. 

The sparse coefficient matrix for an antisymmetric mask is 


np 0 2 y 2 0 0 2 x 2 

A oo 


2 z 

2 y 2 0 0 0 0 

Aoi 


2 zy 

2^ 0 0 2x2y 2 

A 02 


Zzy 2 

Symmetric 2 x 2 0 0 

A io 


2 zx 

2x 2 y 2 0 

A n 


2 zxy 

2 x 4 

A 20 


2 zx 2 


PRECEDING PACT PLANK NOT FILMED 



£i>^s 


This leads to the following reduced matrices. 


np 

2y* 

2x 2 

O 

o 

< 


2z 

2y* 

2y 4 

2x 2 y2 

A«2 

= 

2zy* 

2x? 

2x2^ 

2x 4 

A 2 o 


2zx 2 


2y 2 . A 01 = 2zy 

2X 2 .Aj 0 = 2zx 
S^y 2 . Aj, = 2 zxy 

The evaluation of matrix terms leads to the following values. 


np = (2r + l) (2c + l) = 25 


np np 2 2 2 

I y k = Z \ 2 = Z Z y 2 = 5 £ y 2 

«e=l lc=l x=-2 y=-2 y=-2 


= 50 


np 


k=l 


x k ‘ »k 1 E_ E y* * 5 E 


i,4 = 


170 


k=l 


x= -2 y=-2 


y=-2 


- k y k 2 = Z Z x2 y 2 = 100 


x=-2 y^-2 


Evaluation of the A 00 regression coefficient by Cramer's Rule is 



2 z 

50 

50 


2zy2 

170 

100 

A oo = 

2 zx 2 

100 

170 


25 

50 

50 


50 

170 

100 


50 

100 

170 
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Factoring common terms from rows, expanding the numerator, and 
evaluating the denominator results in 


np / np np \ 

27 2 z - 5 2 zy ! + 2 zx 2 ) I r Norm, 

k=l V k=l k=l / 


where Norm = 175. Recalling that 


np r c 

1=2 I 

k=i x=-r y=-c 


and combining coefficients of like terms , lead to the convolute integers 
A 00 - (- 13z_2_2 + 2 z_ 2 _i + 7z_20 + 2 Z_ 2 i - 13z_ 2 2 


^ 2z_j_2 

N 

1 

-1 

+ 22z_ 10 + 17z_„ + 2 z_ 12 

+ 7z 0 _2 + 

22z 0 _! 

+ 

27z 0 q — 22Zqj 

+ 7z 02 

+ 2Zj_2 + 

17z 1 _ 1 

+ 

22z 10 + 17z 11 

+ 2Z 12 

- 13z 2 _ 2 + 

2*2-1 

+ 

7z 2 q + 2z 21 - 

13z 22 ) t 


Identifying and expressing the C„ coefficients and Norm results in the 
following matrix form for two-dimensional convolute integers for smooth- 
ing, using a 5 by 5 filter on a second order surface. 
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appendix b 





NORMAL. EQUATIONS 
FIRST AND THIRD ORDER SURFACES 
The normal equations for a plane surface constrain m + n ^ 1 


and are 


00 

01 

10 

np 

2y 

Sx 

s y 

2y 2 

Sxy 

Sx 

Sxy 

Sx 2 


For an antisymi etric mask, the normal equations lead to the 
sepa Jble form 


0 0 

Sy 2 0 


Sx 2 I A ni Szx 


where 


A oo = S z/np 


A 10 - 2 zy/L y 2 


A oi = 2 zx/2 x 2 


PREC FDIP'G -• 


r -T FILMED 






- y 

zy 2 

£ xy 

Xy 3 

£ xy* 

X x’y 

Xy* 

£ xy 2 

XxV 

£x 

£xy 

£ x 2 

£ xy* 

£ x*y 

£ x* 

Xxy* 

£ x*y* 

Xx’y 

zy* 

xy* 

£ xy 2 

X y 4 

£ xy* 

£ x*y* 

Zy* 

£ xy* 

XxV 

Exy 

X xy* 

£ x*y 

X xy 5 

£ x*y* 

X x’y 

£xy* 

£ xV 

£xV 

£ x 2 

Z x’y 

£ X s 

£ xV 

X x*y 

X x* 

XxV 

Xx’y 2 

£ x 4 y 


zy* 

£ xy* 

X y 5 

£ xy* 

£ x*y* 

£y* 

£ xy 5 

£ x’y* 

£ xy 2 

£ xy* 

£xV 

£ xy 4 

ZxV 

£ x’y 2 

Xxy 5 

£ xV 

XxV 

£ x’y 

£xV 

X x’y 

£ xV 

£ x’y* 

£ x‘y 

Xx’y* 

£ x’y* 

£ x 4 y* 

£ x* 

£ x’y 

£ X 4 

£ xV 

£ x*y 

£ x 5 

XxV 

£ x 4 y* 

£ x’y 


The antisymmetric property leads to the following sparseness: 


np 0 0 Zy 2 0 Zx 2 0 0 0 0 

0 Zy ! 0 0 0 0 Zy 4 0 ZxV 0 

OOZx 2 0 0 0 0 ZxV 0 Zx 4 


Zy 2 0 
0 0 


0 Zy 4 0 ZxV 0 


0 0 Zx 7 y 2 0 0 0 0 


Zx 2 0 0 ZxV 0 ZX 4 


o n 


0 Zy 4 0 0 

'* 0 ZxV 

Z: 

0 z 


0 Zy* 0 Zx 2 y 4 0 
r '"V 0 ZxV 

















i 



i 


£ 

| The surface order equivalences can now be obtained by considering 

I the orthogonal nature of various groups of columns , which leads to the 

f following separable forms: 


np 

2 y 2 

2x 2 

A oo 


2 z 

Sy2 

2/ 

2xV 

Aq2 

= 

Szy 2 

2x 2 

SxV 

lx 4 

A20 


2 zx 2 

Sy 2 

2y“ 

SxV 

A 01 


Szy 

2/ 

Sy 6 

SxV 

A03 

= 

Szy 3 

2xy 

2 x 2 y 4 

2x 4 y 2 

A 21 


2 zx 2 y 

Sx 2 

2 xV 

2x 4 

a 10 


2 zx 

SxV 

2 x 2 y 4 

SxV 

A 12 

= 

Szxy 2 

Sx 4 

SxV 

2x® 

A30 


2 7.X 6 


SxV . An = Szxy 


Comparing the same partial derivative regression coefficients for 
adjacent surface orders indicates the equivalence to be found in the convolution 
coefficients. The A 00 , A^, and A 20 values are equivalent for second and third 
order surfaces; the A 01 and A 10 terms are equivalent for first and second 
order surfaces. Table 2 represents a general expression for the equivalences 
to be found in the convolution coefficients of adjacent surface orders. 

An alternate definition of the filter coefficients can be expressed as 
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a 

where det(D ) represents the determinant of the cross product matrix of 
order l , 


d*=XS t .XY , 


B 

and B _ .. represents the cofactor of the element in the row, ij column 
of D , as expressed for the second order surface in Appendix A. 

The complementary transpose of an element of a roughing filter can 
be expressed as 


c ji < x >y) ^ Cj. (y,x) = [c* (x,y)f 


From the preceding definition of the convolution coefficients, the above is 
so when 


£ £ 

B o&ji ~ B «j3,iJ ’ 


equality between the cofactors of elements in complementary columns. 

Since the filter masks are symmetric, y is identical with x. Eithe" ignore 
the numerical factors associated with the partial derivatives, Eq. (23), 
or assume equal row column mask spacing, Ar- Ac. Now the normal 
equations indicate that for such masks, the cofactors of elements in comple- 
mentary columns are identical. This leads to the equivalences expressed 
in the lower section of Table 2. 

Furthermore, for these masks 


where either a + i or fi -t j is odd. The convolution coefficients are thus 
concerned with those cofactors where o + i and + j are even. This requires 
the parity of a and i, fi and j to be the same. Now 
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J < ■ ->*W > 4»J i Uttj<W»fc i cyw ,* \»* - -- ’ -jr ; 




rsr* "" %%i* •**•**% 1 


2 ^2°(-l) a x a y^B R .. 

C. (-x,y) tl 

1J det(D) 

Since only parity affects (-1) , this can be 4 \ placed by (-1) ; a and i 
equal parity. This leads to 


Cij(-x.y) = (-l)*C_(x,y) 


Table 3 represents a general expression for the symmetries to be 
found in symmetric filter masks. 
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